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ABSTRACT 

We report a transiting extrasolar giant planet around the star OGLE-TR-10 
(orbital period = 3.1 days), which was uncovered as a candidate by the OGLE 
team in their photometric survey towards the Galactic center (Udalski et al. 
2002a). We observed OGLE-TR-10 spectroscopically over a period of two years 
(2002-2004) using the HIRES instrument with an iodine cell on the Keck I tele- 
scope, and measured small radial velocity variations that are consistent with the 
presence of a planetary companion. This confirms the earlier identification of 
OGLE-TR-lOb by our team and also recently by Bouchy et al. (2004b) as a pos- 
sible planet. Additionally, in this paper we are able to rule out a blend scenario 
as an alternative explanation. From an analysis combining all available radial 
velocity measurements with the OGLE light curve we find that OGLE-TR-lOb 
has a mass of 0.57±0.12 Mjup and a radius of 1.24±0.09 Rjup. These parameters 
bear close resemblance to those of the first known transiting extrasolar planet, 
HD 209458b. 

Subject headings: planetary systems — hne: profiles — stars: evolution — stars: 
individual (OGLE-TR-10) — techniques: radial velocities 

1. Introduction 

Photometric observations combined with radial velocity measurements of stars harbor- 
ing transiting giant planets yield the absolute radii, masses and densities of the planets. 
The ongoing photometric surveys for transiting planets have now delivered large numbers 
of candidates, but those observations alone cannot distinguish between companions that are 
Jupiter-size planets, brown dwarfs, or late type M dwarfs, since they all have similar radii (of 
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order 0.1 Rq). They also cannot distinguish between true transiting planets and so-called 
"blends" — chance alignments or physical triple systems involving an eclipsing binary that 
can mimic planetary transits (Konacki ct al. 2003a; Torres et al. 2004b, 2005). These can 
be very common, particularly in crowded fields. Hence follow-up spectroscopic observations 
are not only necessary to measure the radial velocities and derive the mass of a transiting 
planet, but they are also vital for weeding out astrophysical false positives. 

The first confirmed case of a transiting planet proposed by a photometric survey (Udalski 
et al. 2002a,b) was that of OGLE-TR-56b (Konacki et al. 2003a; Torres et al. 2004a; Bouchy 
et al. 2004b). Three other examples have also come out of the OGLE survey (Udalski 
et al. 2002c, 2003): 0GLE-TR-113b (Bouchy et al. 2004a; Konacki et al. 2004), OGLE- 
TR-132b (Bouchy et al. 2004a; Moutou et al. 2004) and OGLE-TR-lllb (Pont et al. 2004). 
Additionally, the first transiting planet from a wide-field bright-star survey has recently been 
reported (TrES-1; Alonso et al. 2004). This brings the total number to six, including the 
case of HD 209458b, a planet originally discovered in the course of a radial velocity survey 
and only later found to undergo transits (Henry et al. 2000; Charbonneau et al. 2000). 

OGLE-TR-10 was identified as a promising candidate by the OGLE team during their 
2001 campaign in three fields towards the Galactic center (Udalski et al. 2002a). The possible 

planetary nature of its companion based on spectroscopic follow-up was first established by 
Konacki et al. (2003b). In that paper we reported a tentative radial velocity semi-amplitude 
of K = 100 ±43 ms^^, and a mass for the putative planet of Mp = 0.7 ±0.3 Mjup. However, 
the possibility of a blend could not be categorically ruled out at the time due to the small 
number of observations. Recently Bouchy et al. (2004b) also called OGLE-TR-10 a possible 
planet, with K — SI ±25 ms"^ and Mp — 0.66 ± 0.21 Mjup, although these authors were 
still unable to completely exclude a blend scenario because of the insufficient signal-to-noise 
ratio (SNR) of their observations. 

In this paper we present the results from additional spectroscopic monitoring of OGLE- 
TR-10 with Keck I/HIRES on two additional runs, for a total of three seasons (2002-2004). 
We are now able to confirm the planetary nature of OGLE-TR-lOb. Our new radial velocities 
and the parameters of the parent star are presented in Section 2 and Section 3. The OGLE 
light curve sohition and planet parameters are derived in Section 4, where we combine our 
radial velocities with those reported by Bouchy et al. (2004b). In Section 5 we analyze and 
rule out a blend as an alternative explanation. The results are discussed in Section 6. 
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2. Observations 

OGLE-TR-10 was observed spectroscopically with the Keck I telescope using the HIRES 
instrument (Vogt et al. 1994) in July 2002, August 2003, and July 2004. The exposure times 
ranged from 30 to 50 minutes and the wavelength coverage was 3850 A to 6200 A (36 echelle 
orders), at a resolving power of i? ~ 65,000. We collected 10 spectra of which 9 were taken 
with the iodine gas absorption cell (I2) that superimposes a dense forest of absorption lines 
directly on the stellar spectrum in the region from approximately 5000 A to 6200 A (some 
12 echelle orders). The iodine orders had typical SNR of 15 to 25 per pixel. One spectrum 
was taken without the I2 to serve as the template for the iodine velocity reductions. 

For faint stars such as OGLE-TR-10 {V = 15.8, I = 14.9) light losses due to the iodine 
cell typically prevent one from using that technique to determine precise radial velocities 
in the manner done for brighter Doppler targets. However, as described by Konacki et 
al. (2003b) precision as good as 50 ms^^ can still be achieved with the I2 cell on faint 
targets by using a synthetic spectrum instead of an observed spectrum as the template. 
Additionally, in order to monitor changes in the instrumental profile (PSF) that affect the 
velocities significantly, we obtained observations a bright star with the I2 cell on every night 
we observed OGLE-TR-10. Those observations were used to model the PSF and estabhsh 
the parameters that we then applied to OGLE-TR-10. In this way we were able to determine 
radial velocities with typical uncertainties of about 60 ms~^, which are listed in Table 1. 
Our spectroscopic orbital solution is described below in §4. 

3. Pctrameters of the parent star 

The stellar parameters for OGLE-TR-10 were derived from our high- resolution co-added 
Keck spectrum (SNR of 44) with fits to synthetic spectra computed from model atmospheres 
for different compositions based on the ATLAS 9 and ATLAS 12 programs by Kurucz (1995). 
We use a code re- written in Fortran-90 (J. Lester, priv. comm.) and incorporating new 
routines for improved treatment of contributions from various broadening mechanisms, as 
well as updated and expanded opacities and line lists. This code has been tested extensively 
and performs very well for solar-type stars (F-K type) such as OGLE-TR-10. The fits 
between observed and synthetic spectra were made in spectral regions unaffected by the 
I2 lines, and including a large number of metal absorption lines of different ionization and 
excitation states as well as the core and wings of the A4861 H/3 fine. Our effective temperature 
determination, Tgs — 5750 K, has an estimated accuracy of about 100 K (see below). Our 
projected rotational velocity, vsini — 3 kms~^, is good to better than 2 kms~^; we use 
an approach to the fine broadening similar to that described by Fischer & Valenti (2003). 
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Our procedures allow us to derive the metallicity, [Fe/H], and the limb darkening parameter, 
ui — 0.51 ± 0.04, with good confidence. [Fe/H] is estimated to be solar with a 0.2 dex 
uncertainty. The surface gravity is much more difficult to determine reliably. We estimate 
logg = AA^q'q, in which the upper limit is much better constrained than the lower bound, 
as is common in this type of analysis. Given the morphology of evolutionary tracks for 
main sequence stars and the lack of an independent distance estimate for OGLE-TR-10, this 
uncertainty in log g allows us to rule out a giant or subgiant status for the star, but provides 
only a weak constraint on the stellar radius on the main sequence. 

The stellar mass and radius were determined using a stellar evolution code described in 
detail elsewhere (Cody & Sasselov 2002; Sasselov 2003). As evolutionary tracks are nearly 
vertical in the Teflf vs. logg plane for the range of interest for OGLE-TR-10, the largest 
uncertainty in deriving the radius would be an error in T^q. A temperature that is too 
hot would lead to an overestimate of the stellar mass and also a stellar radius that is too 
large (given the stricter upper bound on logg). This in turn would lead to a significant 
overestimate of the derived planet radius {Rp) in the solution of the transit light curve. 
Figure 1 illustrates our efforts to obtain an accurate value of Tcs, and in particular the strong 
sensitivity of Balmer line profiles to the temperature. We point out that LTE radiative 
equilibrium codes like ATLAS are particularly well suited to fit the wings of such lines. 
The cores are usually affected by chromospheres and non-LTE effects, and are expected to 
appear deeper in a solar-type star hke OGLE-TR-10. A fit to the wings of the A4861 H/3 
line requires a careful setting of the contunuum (for better illustration only a fraction of the 
HIRES order is shown in Fig. reffigrspectrum). To guard against problems with contunuum 
setting and line broadening, in deriving Tes we also rely on temperature-sensitive line pairs: 
moderate-strength metal lines having different (often opposite) sensitivity to temperature 
and co-located (within 0.1—0.2 nm) in wavelength (see Gray & Johanson 1991 for the general 
idea and some line pairs). We get T^s = 5800 ± 100 K, in good agreement with the Balmer 
line wings and ncutral-to-ionized metal lines comparisons. 

Thus, based on the overall fit of the good-quality Keck spectrum, OGLE-TR-10 is 
confidently identified as a close analog to our Sun. Our estimates of the stellar parameters 
are significantly different from those reported by Bouchy et al. (2004b) {T^g = 6220 ± 140 
K, log^ = 4.70 ±0.34, [Fe/H] = +0.39 ± 0.14, M = 1.22 ± 0.045 M©, i? = 1.21 ± 0.066 R©). 
However, as can be also seen in Figure 1, our determinations lead to a much better match 
to the overall high SNR spectrum of OGLE-TR-10 from Keck I/HIRES. 



- 5 - 



4. Analysis and results 

Our radial velocities from §2 show clear variations as a function of photometric phase. 
Adopting the orbital period and epoch as reported for OGLE-TR-10 by Udalski et al. (2002a) 
we fitted for a circular orbit solving for the center-of-mass velocity and the semi-amplitude 
K (see Figure 2a). The best-fit value of X = 77 ±23 ms~^ is consistent with the early result 
from Konacki et al. (2003b) and Bouchy et al. (2004b). The rms of our spectroscopic solution 
is 45 ms~^, and from Monte Carlo simulations we find that the probability of obtaining by 
chance a K amplitude as large as we measure is only ~ 4 x 10~^. 

For our final spectroscopic orbital solution we combined our HIRES measurements with 
14 velocity measurements by Bouchy et al. (2004b) ^. The latter were obtained with two 
different configurations of the VLT spectrograph UVES: in the standard sht mode (hereafter 
UVES), and a setup with a fiber link (FLAMES). We considered these as independent data 
sets, and solved for the corresponding velocity offsets relative to HIRES along with the rest 
of the orbital elements. The combined fit gives an improved velocity semi-amplitude of 
= 80 ± 17 ms-\ with an rms residual of 60 ms ^ (see Figure 2b). The corresponding 
false alarm probability from Monte Carlo simulations is ~ 10~^. We use this combined fit in 
the rest of our paper. The best-fit parameters are presented in Table 2, along with our fight 
curve solution based on the OGLE photometry, using the formalism of Mandel & Agol (2002) 
(Figure 3). By combining this information with the stellar parameters described previously 
we derive for the planet in orbit around OGLE-TR-10 an absolute mass of 0.57 ±0.12 Mjup 
and a radius of 1.24±0.09 Rjup. The formal errors include the contribution from uncertainties 
in the mass and radius of the parent star. 



5. False positive rejection 

To test for the possibifity that the radial velocities we measured for OGLE-TR-10 are 
due simply to blending with an eclipsing binary, wc computed the spectral line bisectors 
from our Keck observations as described by Torres et al. (2004a). There is no significant 
variation of the bisector spans as a function of phase (see below), as would be expected if 
lines from another star were causing asymmetries in the profiles of the main star by moving 
back and forth with the photometric period. Similar results were reported by Bouchy et al. 
(2004b). 



^Two additional measurements by these authors were considered by them to be of lower quality and 
rejected. We have done the same here. 
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Next we modeled the OGLE light curve in detail using the procedures described Torres 
et al. (2004b), assuming a configuration consisting of an eclipsing binary blended with the 
main G star forming a hierarchical triple system. Although we were able to achieve an 
excellent fit that is essentially indistinguishable from a true transit light curve, this blend 
model predicts an optical brightness for the primary of the eclipsing binary that is greater 
than the G star itself, which is clearly not observed. We then relaxed the condition that the 
three stars be at the same distance, and considered models with the eclipsing binary in the 
background in order to make it fainter. Here too we were able to find a perfectly acceptable 
fit for an echpsing binary composed of an F9V star and a K5V star located several kpc 
behind the G star (see Figure 4a). The relative brightness of the eclipsing pair in this model 
is only about 4.5% compared to the G star, which would be just below our threshold of 5% 
for detecting lines of another star in the spectra (see Konacki ct al. 2003b). The predicted 
orbital velocity semi-amplitude of the F9 star is 65 kms~^, and its spectral lines should 
show a rotational broadening corresponding to vsini = 20 kms~^ (assuming the spin is 
synchronized with the orbital motion). 

In order to test this blend scenario further we ran extensive numerical simulations 
following Torres et al. (2005) to predict the bisector span and radial velocity variations 
that would be expected from this configuration. In Figure 4 we compare these predictions 
with the observations as a function of orbital phase. The expected bisector span variations in 
Figure 4b are relatively small (less than about 100 ms~^), and are therefore still consistent 
with the measurements, which show no significant variation given typical errors that are 
also about 100 ms~^. Thus, the fact that bisector spans for a transit candidate show no 
appreciable change with phase does not necessarily rule out a blend scenario. However, 
Figure 4c indicates that the expected radial velocity variations are even smaller, and do 
not show the trend with phase displayed by the velocity measurements for OGLE-TR-10. 
As discussed in §4 the latter velocity trend is confirmed and reinforced by independent 
measurements from two different data sets, as reported by Bouchy et al. (2004b). The false 
alarm probability based on Monte Carlo simulations that a set of velocity measurements 
resulting from such a blend scenario would have the phase coherence we observe, with a 
semi-amplitude as large as we observe, is 1.3 x 10~^ from our HIRES measurements alone, 
and 3 x 10~^ using all the RV data. This effectively rules out this blend configuration. 
Experiments show that it is not possible to make the predictions more consistent with all 
of the observational constraints by changing the parameters of the blend model. Additional 
evidence against a blend has been obtained recently by Holman et al. (2005) based on new 
high-quality and high-cadence photometric observations of the star. Their accurate light 
curves in two passbands define the morphology of ingress and egress significantly better 
than the OGLE measurements, and are clearly inconsistent with the light curve from a 
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blend configuration such as we described above, involving a background eclipsing binary. We 
conclude that the observations do not support a blend scenario, and the planetary nature of 
the companion to OGLE-TR-10 is confirmed. 

6. Discussion 

OGLE-TR-lOb is the seventh known extrasolar transiting planet and the fifth to come 
out of the OGLE survey. Admittedly the faintness of the transiting planet candidates from 
the OGLE survey requires some of the largest available telescopes merely to confirm their 
planetary status, and other detailed follow-up studies are very difficult to pursue with current 
instrumentation. Nevertheless, these discoveries around faint stars have been extremely 
important in the field of extrasolar planets. They account for most of the points in the 
current mass-radius diagram for giant planets (see Figure 5), which relates key properties of 
these objects for our theoretical understanding of their structure and evolution. Additionally, 
they have led to the discovery of a new class of "very hot Jupiters" (Bouchy et al. 2004a; 
Konacki et al. 2003a, 2004) with remarkably short orbital periods (1-2 days). The importance 
of faint transit candidates from surveys like OGLE is likely to continue in the near future, 
especially given that the OGLE team has recently released a new set of 40 candidates that 
seem very promising (Udalski et al. 2004). The wide-field surveys, on the other hand, have 
only recently produced their first transiting planet (Alonso et al. 2004) despite having been 
in operation for longer than the OGLE effort. They would provide planets for more detailed 
follow-up studies in the future. 

OGLE-TR-lOb is very similar to HD 209458b (Brown et al. 2001) in terms of its orbital 
and physical parameters. In particular, its orbital period (3.1 days) places it in the "hot 
Jupiter" category (planets with periods of 3-4 days) already well populated from the radial 
velocity surveys. OGLE-TR-lOb has a low average density and hke HD 209458b might 
require an additional heating mechanism to explain it. There are now as many hot Jupiters 
as there are very hot Jupiters among the transiting planets uncovered photometrically, and 
so the apparent lack of longer-period planets in the transit surveys that initially appeared to 
cause some concern (Bouchy et al. 2004a; Pont et al. 2004) no longer seems to be a problem, 
as anticipated by Gaudi, Seager, & Mallen-Ornelas (2004). 

We note, finally, an apparent dichotomy in the mass-radius diagram of Figure 5 in 
that the 4 planets with the longer periods (in the hot Jupiter class) all have small masses 
(~0.7 Mjup), while all the short-period planets (very hot Jupiters) have masses roughly twice 
as large. This trend, noted previously by Mazeh, Zucker, & Pont (2005), now seems to be 
reinforced, and may perhaps be related to issues of survival of planets in very close proximity 
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to their parent stars. 
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Fig. 1. — Portion of the observed (co-added) spectrum of OGLE-TR-10 around the H/3 hne 
with (a) our best-fit synthetic spectrum superimposed (smooth sohd hne) and (b) a synthetic 
spectrum for the stellar parameters derived by Bouchy et al. (2004b). 




Fig. 2. — Radial velocity measurements and fitted velocity curve for OGLE-TR-10, as a 
function of orbital phase. The transit ephemeris is adopted from the photometry (see text). 
Our HIRES velocities are represented with filled circles, and UVES and FLAMES velocities 
by Bouchy et al. (2004b) are shown with open diamonds and open circles, respectively. Panel 
(a) is for the fit to HIRES velocities only (rms of 45 ms~^) and panel (b) is for the fit to all 
available data (rms = 60 ms~^). 
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Fig. 3. — OGLE photometry for OGLE-TR-10 in the / band, with our best fit transit hght 
curve. The resulting parameters are hsted in Table 2. 
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Fig. 4. — Blend model for OGLE-TR-10 in which an eclipsing binary (F9V + K5V) is located 
in the background of the main G star, and has a relative brightness (4.5%) that makes it 
undetectable in our spectra. This particular model assumes a systemic velocity for the binary 
equal to the velocity of the G star, but extensive tests show that this has no effect on the 
conclusions, (a) Fit to the OGLE light curve near the primary eclipse, (b) Measured bisector 
spans as a function of phase, compared to predictions from the blend model. The cusps are 
the result of the large velocity semi-amplitude of the eclipsing binary compared to the width 
of the lines of the G star (see Torres et al. 2005). (c) Our radial velocity measurements shown 
with the predictions from the blend model. The observed velocity trend is not reproduced 
by the model. 
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Fig. 5. — Radii of transiting extrasolar planets plotted against their masses. Jupiter and 
Saturn are included for reference, along with dotted lines of constant density. Data are from 
Brown et al. (2001) for HD 209458b, Torres et al. (2004a) for OGLE-TR-56b, Moutou et 
al. (2004) for OGLE-TR-132b, Pont et al. (2004) for OGLE-TR-lllb, Sozzetti et al. (2004) 
for TrES-1, and this paper for OGLE-TR-lOb. For 0GLE-TR-113b we have combined the 
original observations reported by Bouchy et al. (2004a) and Konacki et al. (2004), and 
derived an improved planetary mass of 1.29 ±0.17 Mjup and a radius of 1.09 ±0.10 Rjup (see 
http : / /www . gps . caltech . edu/~macie j /Planet s/0GLE-TR-113b . html) . 
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Table 1. Radial velocity measurements for 
OGLE-TR-10, in the baryccntric frame. 



HJD 




Velocity 


Error 


(2,400,000+) 


Phase 


(kms^-*^) 


(kms^^) 


52481.7946 


0.706 


-6.130 


0.058 


52483.7578 


0.330 


-6.296 


0.057 


52853.7813 


0.648 


-6.247 


0.048 


52855.8454 


0.314 


-6.340 


0.059 


52864.7786 


0.194 


-6.304 


0.054 


53206.7856 


0.470 


-6.223 


0.097 


53207.7794 


0.790 


-6.152 


0.056 


53208.7772 


0.112 


-6.380 


0.070 


53209.7727 


0.433 


-6.254 


0.056 
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Table 2. Orbital and physical parameters for OGLE-TR-10 and its planet. 



Parameter Value 



Orbital period (days) 3.101386 ± 0.000030 

Transit epoch (HJD-2,400,000) 52070.2223 ± 0.0028 

Center-of-mass velocity (km s"^) -6.250 ± 0.020 

Eccentricity (fixed) 

Velocity semi-amplitude (ms"-*^) 80 ± 17 

Velocity offset between HIRES and FLAMES (ms"^) -14 ± 27 

Velocity offset between HIRES and UVES (ms'^) . . . +218 ± 34 

Inchnation angle (deg) 89.2 ± 2.0 

Stellar mass (Mq) (adopted) 1.00 ± 0.05 

Stellar radius (R©) (adopted) 1.00 ± 0.10 

Fractional radius (ii!pianet/-Rstar) 0.127 ±0.017 

Limb darkening coefficient (/ band) 0.51 ± 0.04 

Planet mass (Mj^p) 0.57 ± 0.12 

Planet radius (Rjup) 1.24 ± 0.09 

Planet density (g cm'^) 0.38 ± 0.10 

Semi-major axis (AU) 0.04162 ± 0.00069 



